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PROFILE-BULK  METHOD  FORMULAS 
FOR  CALCULATING  FLUX  AND  STABILITY 
IN  THE  MARINE  ATMOSPHERIC  SURFACE  LAYER 
AND  A  SURVEY  OF  FIELD  EXPERIMENTS 


Abstract 

A  comprehensive  survey  of  profile  flux  experiments  conducted  in  the  marine  atmospheric  surface 
layer  over  the  last  25  years  is  presented  with  an  extensive  bibliography  and  the  formulas  for  computing 
fluxes  and  stability. 

INTRODUCTION 


A  survey  of  marine  profile  flux  experiments  has  not  appeared  in  the  literature  since  Roll  (1965). 
An  updated  catalog  has  been  prepared  with  an  extensive  bibliography  of  field  experiments,  methods, 
and  instrumentation.  The  constituent  equations  used  to  compute  the  profile  and  bulk  fluxes  are  distri¬ 
buted  over  a  variety  of  sources,  which  contain  different  notations  and  approximations,  and  several  have 
typographical  errors.  A  complete  and  accurate  set  of  formulas  is  presented  for  the  first-time  user. 

SURVEY  OF  MARINE  SURFACE  LAYER  PROFILE  FLUX  DATA 


One  of  the  earliest  experiments  to  be  found  in  the  literature  is  that  of  Wiist  (1920),  which 
reported  wind  profile  measurements  taken  from  a  ship  in  the  Baltic  Sea.  Perhaps  the  earliest  observa¬ 
tions,  as  noted  by  Barenblatt  et  al.  (1975),  were  those  conducted  by  G.  I.  Taylor  in  1913  using  instru¬ 
ments  sent  aloft  with  a  kite  from  the  stern  of  a  whaling  ship  under  sail  off  the  coast  of  Newfoundland. 
A  list  of  marine  surface  layer  profile  measurements  reported  up  to  1962  may  be  found  in  Tables  X  and 
XXIV  of  Roll  (1965).  A  current  survey  of  marine  surface  layer  wind  profile  observations  accompanied 
by  temperature  and/or  humidity  profile  measurements  is  presented  in  Table  1. 


As  indicated  in  the  table,  there  are  two  basic  observational  strategies  for  acquiring  such  profile 
data.  The  first  is  to  employ  a  number  of  identical  sensors  permanently  situated  at  various  altitudes, 
which  acquire  data  on  a  continuous  basis.  The  second  is  to  employ  a  single  sensor,  which  is  sequen¬ 
tially  moved  from  one  altitude  to  another,  acquiring  data  while  stationary  at  the  desired  altitude. 

Our  literature  searches  revealed  the  existence  of  only  about  2,100  h  of  offshore  ocean  profile 
measurements  in  which  all  three  of  the  primary  fluxes  had  been  measured;  85%  of  the  data  had  been 
obtained  in  the  equatorial  region  of  the  Atlantic  Ocean.  It  was  further  determined  that  only  1%  of  that 
data  base  had  been  acquired  under  stable  atmospheric  conditions  and  that  little  data  existed  for  wind 
speeds  in  excess  of  1 2  m/s. 


PROFILE  FLUX  AND  STABILITY  CALCULATIONS 


Encompassed  within  the  profile  method  are  several  semiempirically  derived  profile-flux  relation¬ 
ships  which  have  been  summarized  by  Yaglom  (1977).  The  two  principal  competing  relationship 
schemes  are  those  proposed  by  Dyer  and  Hicks  (1970)  and  by  Businger  et  al.  (1971).  Dyer  (1974)  has 
reviewed  the  differences  the  two  schemes.  Lo  and  McBean  (1978)  found  that  the  two  schemes  could 
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Table  1  —  Marine  Surface  Layer  Wind  Speed  Profile  Observations  Accompanied  by  Temperature 
and/or  Humidity  Profile  Measurements  Reported  Over  the  Last  25  Years 


Source 

Experi¬ 
ment 
(name  or 
location) 

Type  and 

Number  of 

Profile 

Measurement 

Levels 

Altitude 
Cross- 
Section 
Above  Mean 
Sea  Level 
(m) 

Observational 
Strategy 
(M  —  Multiple 
sensors  fixed 
at  measurement 
levels; 

S-Single  sensor 
sequentially 
moved  to 

measure¬ 

ment 

levels) 

Environ¬ 

mental 

Regime 

Measure¬ 

ment 

Plat¬ 

form 

Averag¬ 

ing 

Period 

(min) 

Amount 

of 

Data 

(hours) 

Approxi¬ 
mate 
Range  of 
Wind 
Speeds 
at  10  m 
(m/s) 

Wind 

Speed 

Tempera¬ 

ture 

Humid¬ 

ity 

Deacon  et  al. 
(1956) 

Port  Phillip  Bay 

6 

2 

- 

1.5-13 

M 

Coastal  A 
Inland  Water 

Ship 

30 

49 

1-14 

Takahashi 

(1951) 

Kagoshima  Bay 

5 

5 

5 

0.3-4 

MAS 

Inland  Water 

Boat  A 

Offshore 

Pole 

10 

67 

1-7 

Fleagle  el  al. 
(1958) 

East  Sound 

8 

4 

4 

0.4-4 

MAS 

Inland  Water 

Raft 

60 

30 

3-9 

Bruce  et  al. 
(I960 

Lake  Erie 

6 

6 

- 

0.6-13 

S 

Coastal 

Ship 

30 

20 

3-6 

Deacon 

(1962) 

Port  Phillip  Bay 

3 

2 

- 

3.0-13 

M 

Coastal  A 
Inland  Water 

Ship 

30 

93(?) 

4-15 

Bogorodskiy 

(1964) 

•> 

2 

2 

- 

1.0-10 

M 

7 

Boat 

10 

13 

2-7 

Hoeber 

(1969) 

Equatorial  Atlantic 

4 

4 

4 

1.6-9 

M 

Open  Ocean 

Buoy 

10 

288 

4-10 

Miyake  et  al. 
(1970a) 

Spanish  Banks 

5 

5 

5 

0.6-5 

S 

Coastal 

Offshore 

Mast 

45 

6 

4-9 

Badgley  et  al. 
(1972) 

Indian  Ocean 

6 

6 

6 

1.6—8 

s 

Open  Ocean 

A  Coastal 

Buoy 

40 

79 

2-8 

Oarratt 

(1972) 

Lough  Neagh 

5 

5 

5 

1.0-16 

M 

Inland  Water 

Offshore 

Tower 

10 

<1 

7 

Paulson  et  al. 
(1972) 

BOMEX* 

4 

4 

4 

2.0-11 

S 

Open  Ocean 

Stabilized 

Ship 

48 

113 

2-8 

Hsu 

(1973) 

Fort  Walton  Beach 

6 

3 

- 

0.2-6 

M 

Beach 

Onshore 

Towers 

15 

13 

3-6 

Donelan  et  al. 
(1974) 

Lake  Ontario 

5 

5 

345 

2.4-12 

MAS 

Open  Lake 

Offshore 

Tower 

30" 

132 

1-17 

Dunckel  et  al. 
(1974) 

ATEX" 

7 

4 

4 

1.2-8 

M 

Open  Ocean 

Buoy 

10 

124 

4-11 

Hasse  et  al. 
(1975) 

GATE? 

9 

5 

5 

0.5-8 

M 

Open  Ocean 

Buoy 

10 

6 

2-7 

Peterson 

(1975) 

It  is* 

7 

8 

- 

7.0-123 

M 

Beach  A 

Inland  Coastal 

Onshore 

Tower 

7 

7 

6-9 

Hupfer  et  al. 
(1976) 

Baltic  Sea 

5 

5 

- 

1.5-7 

S 

Beach 

Onshore 

Pole 

30 

4 

-7 

KrOgermeyer 

(1976) 

ATEXb 

7 

4 

4 

1.2-8 

M 

Open  Ocean 

Buoy 

10 

124 

4-11 

Vugts  A 
Businger  (1977) 

Schiermonnikoog 

Island 

5 

5 

- 

0.5-13 

M 

Beach 

Onshore 

Towers 

60 

45 

4-8 

Hasse  et  al. 
(1978b) 

GATE' 

7 

5 

5 

1.1-8 

M 

Open  Ocean 

Buoy 

10 

1.419 

1-12 

Thomson 

(1979) 

Ris* 

4-1 

2-3 

1.0-27 

M 

Beach  A 

Inland  Coastal 

On-  A 

Off-shore 

Towers 

10 

32 

2-7 

Blanc 

(1981) 

San  Nicolas  Island 

2 

2 

2 

9.2-18 

... 

M 

_ 

Coastal 

Onshore 

Tourer 

30 

136 

2-17 

••OMEX  —  Barbados  Oceanographic  and  Meieorotoflcal  Experiment,  1969. 

"ATEX  m  Atlantic  Trad*  Wind  Experiment,  1969. 

cOATE  m  Global  Atmospheric  Research  Program  Atlantic  Tropical  Experiment,  1974. 

11  JO- mm  turning  average  updated  ever  10  min.  2 
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yield  differences  as  large  at  40%  in  the  estimated  fluxes;  they  suggested  that  the  discrepancies  could  be 
resolved  by  setting  the  von  Karman  constant  equal  to  0.40,  instead  of  the  recommended  0.35,  when 
using  the  Businger  et  at.  scheme.  This,  in  effect,  renders  the  Businger  et  al.  scheme  equivalent  to  that 
of  Dyer  and  Hicks.  A  justification  for  this  approach  was  subsequently  provided  by  Wieringa  (1980). 
For  additional  information  see  Wyngaard  et  al.  (1982),  Wieringa  (1982),  and  Gill  (1982).  An  interna¬ 
tional  turbulence  measurement  comparison  conducted  over  land  and  reported  by  Francey  and  Garratt 
(1981)  proposed  a  von  Karman  constant  of  0.38  (±0.04).  An  alternative  analysis  of  the  same  compari¬ 
son  presented  by  Dyer  and  Bradley  (1982)  suggested  a  constant  of  0.40.  A  rationale  for  the  minimum 
height  of  the  lowest  measurement  level  and  the  separation  between  adjacent  levels  is  presented  in  Blanc 
(1983a). 

At  altitudes  _z  —  6,  12,  and  24  m,  for  example,  30-min-averaged  measurements  are  made  of  the 
air  temperature  ( T„ )  in  °C,  the  dew-point  temperature  ( Td„ )  in  °C,  and  the  wind  speed  ( u„ )  in  m/s. 
The  30-min-averaged  barometric  pressure  ( P )  in  pascals  is  measured  at  the  middle  altitude  and  is  used 
for  all  levels.  For  simplicity,  the  overbar  notation  used  to  indicate  time-averaged  measurements  will  be 
implied  but  not  indicated.  Thus,  at  z  —  24  m  (zu)‘  Tm,  7^24*  «24i  at  z  —  12  m  (z^):  Ti2,  Td]2,  ui2, 
P,  and  at  z  —  6  m  (z6):  f6,  Td 6,  m6.  If  the  dew-point  temperature  is  not  measured  at  each  altitude, 
then  the  wet-bulb  temperature  or  relative  humidity  should  be  measured  instead. 

The  saturated  vapor  pressure  ( esn )  in  pascals  at  each  altitude  (z„)  is  calculated  by 


ejP.T) 


p.Q  Q  §+ Q  ^ 


where 


373.16 
T„  +  273.16' 


a2  -  at  -  1, 

a3  -  1  _  — 

«4  -  (I0'2*>)  -  1, 
as  -  (10a3*J)  -  1, 

the  Goff-Gratch  humidity  formulation  constants  are 

6,  -  -3.49149,  64  “  -7.90298, 

62  -  11  344,  65  -  8.1328  x  irJ, 

63  -  5.02808,  b6  -  -1.3816  x  10"7; 

and  n  —  6,  12,  and  24. 


The  water  vapor  pressure  (e„)  in  pascals  at  each  altitude  (z„)  is  calculated  by 

e„(P,  Td)  -  Pa\l  x  iofl^+fl<*5+»s\ 


where 


373.16 

Ta.  +  273.16’ 
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a 2  through  o,  are  calculated  in  the  same  manner  as  above,  and  n  -  6,  12,  and  24.  In  the  event  that 
the  wet-bulb  temperature  (7^,,)  in  °C  or  the  relative  humidity  (RH„)  in  %  is  measured  at  each  altitude 
instead  of  the  dew-point  temperature  ( Td„ ),  then  the  water  vapor  pressure  can  be  calculated  by 

eJP.T.e,,^)  -  esn  -  6.6  x  I0“4[1  +  (1.15  x  10"3  Twbn))P(T„  -  Twbn) 
or 


c„(c„RH) 


100  ’ 


The  specific  humidity  (?,) 


in  kg/kg  at  each  altitude  (z„)  is  calculated  by 


4*(P.e) 


0.622  eK 


P  -  (0.378  e„)  ‘ 


The  absolute  humidity  or  water-vapor  density  (pv„)  in  kg/m3  at  each  altitude  (z„)  is  calculated  by 

pv„(T,e)  -  c«[  -I-  273.16)  ’ 

where  Ry  is  the  water-vapor  gas  constant,  4.615  x  102  J/(kg  K).  The  mixing  ratio  (rj  in  kg/kg  at  each 
altitude  (z„)  is  calculated  by 


0.622  e„ 
P~en 


The  relative  humidity  (RH„)  in  percent  at  each  altitude  (z„)  is  calculated  by 


RH„(e,eJ)  -  —  x  100. 


In  each  case  n  -  6,  12,  and  24. 


The  virtual  temperature  (Tv„)  in  °C  at  each  altitude  (z„)  is  calculated  by 
Ty„(T,q)  -  {(r„  +  273.16)  x  [1  +  (0.608?,)])  -  273.16. 

The  potential  temperature  (0„)  in  °C  at  each  altitude  (z„)  is  calculated  by 

9„(T,z)  -  T„  +  (0.0098z„). 

The  virtual  potential  temperature  (0V„)  in  °C  at  each  altitude  (z„)  is  calculated  by 

9vn(Tv,z)  -  Ty„  +  (0.0098z„). 

Again,  n  -  6,  12,  and  24. 

The  geometric  mean  height  (gmh)  is  the  geometric  middle  altitude  when  vertical  distance  is 
represented  on  a  log  scale.  Because  the  three  measurement  altitudes  (z6,  z)2,  z24)  were  selected  such 
that  (In  z24  —  In  Z|2)  ■»  (In  z)2  -  In  z6),  zi2  is,  also,  the  gmh: 

gmh  -  (z6z 24)l/2  -  Z|2. 

Since  profiles  tend  to  vary  approximately  logarithmically,  the  gmh  of  the  profile  measurement  levels  is 
the  altitude  for  which  the  profile  measurements  are  typically  most  valid. 

The  Richardson  number  stability  (Rigmh),  a  dimensionless  quantity,  valid  for  the  geometric  mean 
height  of  the  component  measurement  altitudes  is  calculated  by 
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^*gmh  “  ^*12 


where 


(Tm  +  273.16)  x 


7v24  ®v6 


dz  gmh  ln(z24/z6)  ’ 

di/  __  “24  -  Uf, 

dz  ~~  gmh  ln(z24/zg)  ’ 

gmh  —  z  i2, 

and  g  is  the  acceleration  constant  due  to  gravity,  9.8  m/s2.  The  approximations  for  the  partial  deriva¬ 
tives  are  accurate  to  about  ±  3%. 

If  Ri|2  <  0,  the  atmosphere  is  unstable  and  the  Monin-Obukhov  stability  (£|2)  at  z  —  12  m  Czi2) 
is  calculated  by 

£,2(-0.01  <  Ri  <  0)  sb  1.3  Ril2, 


where 


£12(—  1.5  <  Ri  <  -0.01)  =  -10f| 

c,  -  2.8443  x  10“2, 

c2  -  0.96125[log(— Ri|2)l, 

c3-  1.3655  x  ir3llog(-Ril2)l2. 


{l2(-oo  <  Ri  <  -  1.5)  =  1.05  Ri,2. 

The  middle-range  (l2  approximation  is  a  least-squares  polynomial  fit  to  the  reverse  solution  of  the  Bus- 
inger  et  al.  (1971)  equation  for  Ri  ({  <  0).  All  three  of  the  approximations  are  accurate  to  about 
±2%.  The  quantity  £  is  of  questionable  validity  for  Ri  <  -2. 

The  Monin-Obukhov  length  (L)  in  meters  is  then  calculated  by 

L  -  -’2- 

£.2' 

Once  L  is  calculated  for  one  altitude,  the  Monin-Obukhov  stability  (£„)  for  any  altitude  (z„)  can  be 
computed  by 

«■(*>  S  f 

When  Ri|2  <  0,  then  {i2  <  0.  Under  near-neutral  conditions  ({|2  =  0)  the  various  profiles  exhibit  a 
linear  form  when  plotted  on  a  semilogarthmic  plot.  Under  non-neutral  conditions  the  profiles  need  to 
be  corrected  for  stability-dependent  curvature  (see  Fig.  1)  to  properly  compute  the  fluxes.  The 
stability-  and  altitude-dependent  profile  corrections  for  potential  temperature  for  specific  humidity 
ty4),  and  for  wind  speed  («#»„)  are  calculated  by 

+•*  ({  <  0)  -  (£  <  0)  -  In  -^^-1 


BaKSjygjaaa 


3CE3aE3Kze?3BS^^ :>y:--i--..:v'v'.v.-:'.q 
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Fig.  I  —  Typical  stability-dependent  curvatures  in 
a  semilogarithmic  profile  representation  produced 
by  non-neutral  conditions.  The  scaling  parame¬ 
ter,  in  this  case  the  friction  velocity  (u,),  is  a 
convenient  shorthand  for  expressing  the  profile 
and  is  inversely  proportional  to  the  stability- 
corrected  profile  slope.  The  roughness  length 
(z0)  is  defined  as  the  altitude  at  which  the  wind 
speed  would  be  zero  if  the  stability-corrected 
wind  profile  were  projected  downward. 


and 


where 


({  <  0)  -  2  In 


I  +  x. 


+  In 


I  +  (*„)* 


-  [2  arctan  x„\  +  y , 


arctan  x„  is  expressed  in  radians,  and 


*  -  Cl  -  15{ .)* 

y.  -  (1  - 

n  —  6,  12,  and  24. 


If  Ri|2  «■  0,  the  stability  is  said  to  be  neutral,  {|2  “  0,  and 

-  0)  -  -  0)  -  -  0)  -  0. 

If  Ri(2  >  0,  the  atmosphere  is  stable  and  the  Monin-Obukhov  stability  (Cn)  at  z  -  12  m  (zt2)  is 
calculated  by 

,  „  ~di~  <<*22  -4<M3)i/2 

Ci2(0.2  >  Ri  >  0) - - - \r. - — — , 

zdi 

where  d\  through  </j  are  the  components  of  a  quadratic  reverse  solution  of  the  Businger  et  al.  (1971) 
equation  for  Ri({  >  0),  in  which 

di  -  (22.09  Riu)  —  4.7, 
d2  -  (9.4  Ri,j)  -  0.74, 


dj  “  R«i2- 


V.V  V  *.•/.  .  \ 

.  *  ,*  ,.»4  •  ■  *.•*•*. 


NRL  REPORT  8647 

The  quantity  (  is  undefined  for  Ri  >  —0.2.  The  Monin-Obukhov  length  (L)  in  meters  is  calculated, 
as  before,  by 

l-4*-. 


Again,  once  L  is  calculated  for  one  altitude,  the  Monin-Obukhov  stability  ((„)  for  any  altitude  (z„)  can 
be  computed  by 

l.M  S  f . 

When  Ri|2  >  0,  then  {12  >  0.  The  stability-  and  altitude-dependent  profile  correction  for  potential 
temperature  for  specific  humidity  tyr,),  and  for  wind  speed  (</>„)  are  calculated  by 

>  0)  -  *,„((  >  0)  -  Z^JL 


where 


*«,({  >  0)  -  -4.7£„, 


n  —  6,  12,  and  24. 


The  mean  stability-corrected  slope  of  the  potential  temperature  profile  <St>  is  calculated  by 

^  ^  £»(6,I2>  +  £*(12.24)  +  2  £#(6. 24) 


where 


£»(».*) — 


On  2^  -  +tm)  -  On  r,  -  *,„) 


and  n,m  —  6,12;  12,24;  and  6,24.  The  mean  stability-corrected  slope  of  the  specific-humidity  profile 
<£,>  and  the  mean  stability-corrected  slope  of  the  wind-speed  profile  <S„>  are  calculated  by  insert¬ 
ing,  respectively,  q  for  all  0’s  and  u  for  all  0’s  in  the  above  two  equations.  The  profile  slopes  are  util¬ 
ized  here,  instead  of  the  more  conventional  profile  gradients,  to  eliminate  the  confusion  that  exists  in 
the  literature  between  the  contradictory  Western  and  USSR  definitions  of  gradient  and  lapse  rate. 

The  scaling  potential  temperature  (0  *)  in  *C,  the  scaling  specific  humidity  (q«)  in  kg/kg,  and  scal¬ 
ing  wind  speed  («*)  in  m/s  (see  Fig.  1)  are  calculated  by 

a  _  * 

*  0.74  <S#>  ’ 

q*~  0  74  <Sf  >  ' 


<S„>  ’ 

where  k  is  the  von  Kirmdn  constant  of  0.40.  The  scaling  wind  soeed  parameter  (u„)  is  more  com¬ 
monly  called  the  friction  velocity.  The  roughness  If  >th  (zq)  in  iters  is  computed  from  the  wind- 
speed  profile  by 

(In  *ti-*,ij)- 1 

*0“  * 


•  s'.  --'S' 
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The  density  of  moist  air  (p)  in  kg/m3  for  all  altitudes  is  computed  by 


p(P,Tv ) 


(3.4838  x  10“3)/> 


rv  ,2  +  273.16  • 

The  specific  heat  of  moist  air  at  constant  pressure  ( Cp )  in  J/  (kg  K)  for  all  altitudes  is  computed  by 

Cp(i j)  =  1.004(1  +  0.9(?12)]  x  103. 

The  latent  heat  of  water  vaporization  (Z,v)  in  J/kg  for  all  altitudes  is  computed  by 

LV(D  =  4.1868(597.31  -  0.56525r12)  x  103. 


Employing  a  negative  flux  value  to  indicate  a  downward  direction  and  a  positive  flux  value  to  indi¬ 
cate  an  upward  direction,  the  sensible  heat  flux  (Hs)  in  W/m2  is  computed  by 

Hs  -  - pC„0  ,u„ 


The  water- vapor  flux  (£)  in  kg/(s  m2)  is  computed  by 

E  -  -p  q*u,. 

Frequently  the  moisture  flux  is  presented  in  terms  of  the  latent  heat  flux  ( HL )  in  W/m2, 

Hl  -  E  Iv 

The  momentum  flux  (M)  in  N/m2  is  computed  by 

M  -  -pu„2. 

Frequently  the  mechanical  flux  is  represented  by  the  stress  (r)  in  N/m2, 

t5-M 


BULK  FLUX  AND  STABILITY  CALCULATIONS 

Encompassed  within  the  bulk  method  are  more  than  20  different  semiempirical  bulk-flux  relation¬ 
ship  schemes.  Partial  summaries  may  be  found  in  Pond  et  al.  (1974),  in  Friehe  and  Schmitt  (1976), 
and  in  Kondo  (1977).  The  various  types  of  proposed  schemes  run  the  complete  gamut  in  terms  of 
complexity  and  sophistication.  For  example,  where  Friehe  and  Schmitt  use  single  statistical 
conglomerate  coefficients,  Kriigermeyer  (1976)  employs  stability-dependent  coefficients,  and  Liu  et  al. 
(1979)  utilize  a  computational  iteration  technique  requiring  five  input  parameters  for  twelve  dependent 
variables  in  five  equations  which  reduce  to  three  simultaneous  equations  of  three  unknowns.  At  the 
present  time  there  is  no  single  universally  accepted  bulk-flux  scheme.  Depending  upon  which  scheme 
is  used,  the  same  input  data  can  yield  estimated  fluxes  differing  by  as  much  as  100%. 

Admittedly,  because  of  the  wide  variety  from  which  to  choose,  the  selection  of  a  bulk  scheme  to 
accompany  the  profile  method  is  of  necessity  somewhat  arbitrary.  Here  1  employ  the  sensible  heat  flux 
and  latent  heat  flux  schemes  proposed  by  Friehe  and  Schmitt  (1976),  as  used  in  Friehe  and  Pazan 
(1978),  and  the  momentum  flux  scheme  proposed  by  Smith  (1980a).  1  have  selected  these  schemes 
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because  of  their  simplicity  and  statistically  broad  data  bases.  Other  types  and  forms  of  profile-bulk 
combinations  have  been  proposed  by  Okamoto  et  al.  (1968),  Fujita  (1978),  Itier  (1980),  and  Wang 
(1981). 

In  addition  to  those  measurements  made  for  the  profile  calculations,  the  bulk  calculations  require 
a  water-temperature  measurement  ( Tw)  in  °C,  made  at  or  just  below  the  sea  surface.  The  potential 
temperature  (01O)  in  °C,  the  absolute  humidity  or  water-vapor  density  (pvl0)  in  kg/m3,  and  the  wind 
speed  («io)  in  m/s  at  the  standard  bulk  measurement  altitude  of  z  —  10  m  (z10)  can  be  computed  from 
the  profile  measurements  at  z  -  6  m  (z6)  and  z  -  12  m  (z]2)  by 

0io  ~  0]2  —  0.263  (0i2  —  6^), 

PvlO  ~  Pvl2  ~  0.263  (pvi2  —  Pv6^ > 


“to  ~  «12  -  0-263  («,2  -  u6), 


where 


In  Z12  -  In  zip 
In  Z]2  -  In  z6 


=  0.263. 


The  mean  vertical  velocity  (w) -temperature  (0)  covariance  <w'0'>  in  (m  K)/s  is  computed  by 

<<u'0'>  ~  0.002  +  [ChU\0(Tw  ~  0jo)J> 

where  CH  is  the  bulk  sensible  heat  transfer  coefficient  at  the  10-m  altitude,  0.92  x  10-3. 

The  sensible  heat  flux  (Hs)  in  W/m2  is  computed  by 

Hs  -  pCp<(o'0'> , 

where  the  density  of  moist  air  (p),  the  specific  heat  of  moist  air  at  constant  pressure  (Cp),  and  the  flux 
direction  sign  convention  are  the  same  as  used  in  the  profile  calculations. 

The  water-vapor  pressure  near  the  water  surface  iew)  in  pascals  is  calculated  by 


where 


ew(P,Td)  -  Pa?3  x 


373.16 


1  Tdw  +  273.16’ 

the  barometeric  pressure  ( P )  is  in  pascals,  a2  through  a5  and  b\  through  />6  are  the  same  as  for  the 
profile  calculations  of  vapor  pressure,  and  the  dew  point  near  the  water  surface  ( Tdw )  in  °C  is  estimated 
by 

~  Tw, 

assuming  the  relative  humidity  at  or  near  the  sea  surface  is  100%.  The  absolute  humidity,  or  water- 
vapor  density,  at  the  sea  surface  (pvw)  in  kg/m3  is  then  calculated  by 

pvw(T,e)  -  R  {Tw  +  211.\b)  ’ 
where  Rv  is  the  water- vapor  gas  constant,  4.615  x  102  J/(kg  K). 
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The  mean  vertical  velocity  (tu)-absolute  humidity  (pv)  covariance  <at'p's>  in  (m  kg)/(s  m3)  is 
calculated  by 

<«u'pv>  =*  CgUl0(pvw  —  pvio). 

where  CE  is  the  bulk  water-vapor  transfer  coefficient  at  the  10  m-altitude,  1.32  x  10~3.  The  water- 
vapor  flux  (£)  in  kg/fs  m2)  is  computed  by 

E  ”  <o»'pv>  • 

As  before,  the  latent  heat  flux  (HL)  in  W/m2  is  computed  by 

Hl  -  E  Ly, 

where  Z,v  is  the  same  latent  heat  of  water  vaporization  employed  in  the  profile  calculations. 

The  momentum  flux  (A/)  in  N/m2  is  computed  by 

M  “  ~pCpUio, 

where  CD  is  the  bulk  aerodynamic  drag  coefficient  at  the  10-m  altitude, 

CD  -  (0.61  +  0.063  u10)  x  1<T3. 

As  before,  the  stress  (r)  in  N/m2  is  computed  by 

r  =  —  M. 


The  scaling  potential  temperature  (9  ,)  in  °C,  the  scaling  specific  humidity  (qj  in  kg/kg,  and  the 
scaling  wind  speed  or  friction  velocity  (u.)  in  m/s  are  calculated  by 

pu*Cp 

~E 


The  Monin-Obukhov  length  ( L )  in  meters  is  calculated  by 

~(9l0  +  273.1 6)  u  i 
L  ”  gk<m'9’>  ’ 

where  the  acceleration  due  to  gravity  (g),  9.8  m/s2,  and  the  von  Karman  constant  ( k ),  0.40,  are  the 
same  as  used  in  the  profile  calculations. 

The  Monin-Obukhov  stability  ((|2)  for  the  geometric  mean  height  of  the  profile  measurements  at 
z  -  12  m  (gmh  —  Z|2)  is  computed  by 


If  C,2  <  0,  the  unstable  condition,  the  Richardson  number  stability  (Ri|2)  is  computed  by 


Rii2({  <  0) 


0.74  {12 


(1-15{,2)i/2 
(I  -  9<I2)1'2  ' 
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If{i2-  0,  the  neutral  condition,  the  Richardson  number  stability  (Ri!2)  is  also  zero. 


If  £12  >  0,  the  stable  condition,  the  Richardson  number  stability  (Ri^)  is  computed  by 

(0.74  +  4.7{,i){,! 

r,i!<{  >  0)  -  -r+i.n,:)1- ' 


The  roughness  length  (z0)  in  meters  is  computed  from  *fiu |2  (via  £t2)  and  ut  in  the  same  manner 
as  for  the  profile  method. 
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